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ABSTRACT 

A mathematical model is constructed by using the thermal properties and thermal mathematical equations. For 
a particular geometry of a compound parabolic concentrator such as acceptance angle 30°, concentration ratio 2, length 
1.6 meter, the thermal efficiency and performance is calculated by using the mathematical model. The solar radiation 
fall on the compound parabolic concentrator is also calculated. The thermal performance, efficiency and solar radiation 
is calculated for Patan, Gujarat, India (23.4°N, 72°E). The effects of air mass flow rate, Wind Velocity and solar 
radiation on the thermal performance of compound parabolic concentrator is calculated theoretically which is close 
agreement with experimental data that indicate thermal efficiency increase as mass flow rate increase and gain in 
temperature decreases. As wind velocity increases, the overall loss coefficient increases so that thermal efficiency 
decreases. 
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INTRODUCTION 

Compound parabolic concentrator is non-focusing type concentrating collector and provides moderately 
solar radiation concentration although less than a focusing collector it does not require continuous tracking of the 
sun but it need necessities only a few tilt adjustment per year. The concentration ratio (C) is the ratio of area of 
aperture to the area of receiver, which is generally ranging from 2 to 9. The acceptance angle is denoted by 0 a and 

calculated as € = — - As having higher acceptance angle the efficiency of accepting diffuse radiation is much 

larger than the conventional concentrators. Compound parabolic concentrator consists of two parabolic reflector, 
glass cover, flat plate absorber as shown in figure [lj.Rays in the central region of the aperture reach the absorber 
directly where as those near edges undergo one or more reflection before reaching the absorber. Reflector is 
selected in such a way that it reflects the ray at high rate. Reflector should not consist higher absorptivity. In 
opposite receiver is selected in such a way that it have higher absorptivity and lower reflectivity. It also consists 
higher thermal conductivity so that it transfers heat to the fluid at high rate. Most receiver material does not have 
high absorptivity so they need to be covered with special solar selective surface coatings. The rays inside the CPC 
have lower wave length so there are possibilities of loss of these infrared rays. The glass cover does not allow to 
these rays to go outside, so it reduces the loss and also protects CPC from dust and other particles. The properties 
of air are shown in table 1. The characteristics of CPC are shown in table 2[15]. The schematic fig of compound 
parabolic concentrator is shown in figure 1. 
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REVIEW OF LITERATURE 


A solar air heater which heats the air by solar energy it have many application in drying agricultural products such 
as seeds [1], fruits and vegetables [2-3]. Solar air-heaters are utilised for heating buildings with auxiliary heaters to save 
energy in winter-time [4]. Garg et al. [5] concluded that the multiple-pass air heater performs well for low flow rates and 
large plate lengths. 

Wijeysundera el al. [6] perform experiments and concluded that in two pass mode (inlet air flows first above the 
receiver plate then under it.) increase efficiency up to 10-15%. Ezeike [7] concluded that the a triple pass Hate collector 
gives the average efficiency of 70-80% with air velocity up to 3.5 m/s. Yeh el al. [8] and El-sebili et al. [9] designed double 
flow solar air heater with fins attached over and under the absorbing plate. Tchinda and Ngos [10] developed equation for 
calculating thermal performance of CPC collector with a flat side receiver, with various dimensions and concluded that 
efficiency increase as mass flow rate increase and gain in temperature decreases. 

METHODOLOGY 

Calculation of Solar Radiation 

Solar radiation can be found by calculating the beam radiation and diffuse radiation. These radiations can found 

by [11]. 

Calculation of Diffuse Radiation 


The diffuse radiation can be found by [13] 




Where Pi/Po is the pressure ratio at location concentrated relative to the standard atmospheric pressure given as 

p , 

-^= exp C— 0.0001134 


Where H ll]t is altitude in (m) from sea surface. 

The extra-terrestrial solar intensity Al, the atmospheric extinction coefficient B, and sky-diffused factor. Cl, were 
calculated by Joudi (1988) [12] for any day by the following equations: 

WD 

= 1153 X (1 +- 0.066* ena(360 X — J) 

3 0 


B = 0.175 X [(1 — 0,2 X cos(0,93 X NO)] - 0 ,0045 X [1 -cos(l.S6 X 1YU}] 
= 0,0965 x [l - 0.42 x oos(^- x W)] - 0.0075 x [l - ccs(1.95 x ND}] 


Where, ND=number of day in year. 

The sun altitude angle a is given by (Lunde, 1980) [13] as: 
a = sin - L (cos © x cos £ X cos vr 4- sin £ X sin 0) 
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The declination angle 
36U ,■ 

6 = 23.4-5 sin(— x (2 35 + ND}) 

365 

Where, 

* = latitude angle (for Patan * = 23.8667) 
w = the hour angle, this angle is expressed in degree as: 
w = 15 X (solar time — 12) 

Where the solar time calculated as: 

Solar time= standard time + (E/60) - ((L st -L local )/15) 

Where, 

E = correction time which is in hour, which is taken as 0.020833 
L st = standard meridian for local time zone (for india it is 82.5°) 

Liocai= longitude of the location in degree (for patan it is 72.0167°) 

Calculation of Beam Radiation 

The beam radiation is given by 

1 3 = [ptf X 

Where, cos is calculated by [14] : 

cos A = [(sinfi x sin. 0 x cos/ 2 ) - (sitiiS x cos 0 x sin/? x cos 7 ) + (cos £ x cos/? x cos 0 x cosm) + 

(cos 6 X sing? X sin/? X cos y X cos m) + (cos £ X sin/? X sin . y X sin w)] 

Where 

r= surface azimuth angle, it is on horizontal plane between the line due south and the horizontal projection on the 
normal to the inclined plane surface. 

P = slop of the collector (23° for summer and 46° for winter) 

The diffuse radiation on tilted collector is given (Lunde, 1980) [13] as: 

[ 1 — cos£? I 
(A. X (1 + eos/?)/2)) + (S X (Cl + sines) X ( — -)) 

Where 

S= ground reflectivity, which is 0.2 
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Total Solar Radiation 

The total solar radiation is calculated by 

It = "l" J b 


Calculation of Overall Loss Coefficient 

Compound parabolic concentrator consists various heat losses. The wind velocity also affects the performance of 
the CPC. If the velocity is higher, then the loss is increased and so that lower efficiency is obtained. 

Calculation of Wind Heat Transfer Coefficient (h„): 

The convective heat loss coefficient can calculate by Duffie and Beckman [14] as: 


Reynold number = 




Where, 

p a ; f = Density of air, kg/m 3 
(t er = Dynamic viscosity of air. Ns/m 2 


velocity of wind, m/s 


K % w a +h 


Where, w a =width of aperture 
h =height of compound parabolic concentrator 
Nusselt number: 

Nu= 0.40+0.54Re° 52 for 0.1<Re<1000 
=0.30xRe° 6 for 1000<Re<50000 

Nussslt number = x 

Where, 

Ji, v ,= wind loss 

fc s ,'r = thermal conductivity of air 

The heat loss is calculated in two parts. The heat loss and overall loss coefficient are calculated by Duffie and 
Beckman [14]: 

Heat loss From Absorber Plate to Cover Plate 

• Conduction loss 
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• Radiation loss 

Heat Loss from Cover Plate to Sky 

• Convection loss due to wind 

• Radiation loss 

Conductive Heat Loss from the Receiver Surface to Cover Plate 

Qcoprtf = Kair x 4 X (T s - T c ') 


Radiation Heat Loss from Receiver Surface to Cover Plate 


■ - 


4 Xflxt Tg-T*) 



Convective Heat Loss from Cover Plate to Air 

Qctufv- = 4 x 4 x fJc ~ 


Radiation Heat Loss from Cover Plate to Sky 

Qrad: ■ = 4 x a x (Tjf - 

Overall Loss 

= QcONIf ^QsADi. "f" QcOIW ^QuADI 
Overall Loss Coefficient 


4 = 


QlDS5 

4 x (T r - t s 3 


Thermal Efficiency 

The thermal efficiency depends on various factors such as overall loss coefficient (U L ), solar radiation (I), aperture 
area (A a ), receiver area (A R ). The thermal efficiency is calculated by Duffie and Beckman [14] as: 

Qu = * X 4 x (H - £ x U L x (Tj, - r A )> 

Where, 

H = solar radiation x reflectivity of reflector x absorptivity of absorber 
4= Area of aperture without shading 

4 ,T a ,T r = Inlet, Ambient, Receiver temperature respectively 

4 = F' x F 
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F = — 

u r 


m' x C p r 


Rise in Temperature 

Qu 


-AT = ■ 


. X C- 


Outlet Temperature 

T b = T } +- NT 


1 

™- x cy J 


Efficiency 


1 = 


Qv 

1 X A a 


RESULTS AND DISCUSSION 

Figure 2 shows Time vs. solar radiation, it shows that solar radiation is change throughout the year. It is maximum 
in March to June and minimum in September to December. Fig 3 shows wind velocity vs. overall loss coefficient, it is 
concluded that as wind velocity is higher, then loss of heat is at high rate, which reduces the thermal efficiency. Figure 4 
Shows Day time vs. Temp, it is concluded that outlet temp is maximum at noon and minimum at morning. Figure 5 shows 
the outlet temperature vs length for different mass flow rate.it shows that the outlet temperature is higher as length 
increases, it also shows that at mass flow rate m=0.012 kg/s higher outlet temperature is obtained. 

Figure 6 shows mass flow vs outlet temperature for different length L= 1 .6, 2, 2. 2 meter. It shows as mass flow rate 
increase the outlet temprature decreases. Figure 7 shows mass flow vs. efficiency for different velocity lm/s and 8m/s. it 
shows that due to higher velocity heat loss is higher, so that at higher velocity the efficiency is lower. And thermal 
efficiency increases as mass flow increases. Figure 8 shows Day time vs. Efficiency, it is concluded that efficiency is also 
varies with time, which is higher at noon and lower at morning. 

CONCLUSIONS 

As we increase the mass flow rate, outlet temperature decreases and thermal efficiency increases. As wind 
velocity increases, the overall heat loss increases, which reduces rise in outlet temperature so that efficiency also decreases. 
This method is useful to calculating the effect of mass flow rate, wind velocity on thermal efficiency of Compound 
parabolic concentrator. 
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Figure 1: Schematic View of Compound Parabolic Concentrator 




Figure 3: Wind Velocity vs. Overall Loss Coefficient 
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OUTLET TEMPERATURE VS LENGTH 
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Figure 5: Outlet Temperature vs Length (for Mass Flow Rate m=0.02, 0.016,0.012 kg/s) 


MASS FLOW VS OUTLET TEMPERATURE 
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Figure 6: Mass Flow vs Outlet Temp.(for Length L= 1.6, 2, 2.2 meter) 
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Table 1: Properties of Air 


Parameter 

Unit 

Value 

Density 

Kg/m” 

1.127 

Dynamic viscosity 

Ns/m 2 

0.00001983 

Thermal conductivity 

W/mk 

0.0271 


Table 2: Characteristics of Compound Parabolic Collector 


Parameter 

Unit 

Value 

Cover absorptance 

- 

0.05 

flat plate absorber absorptance 

- 

0.95 

Cover transmittance 

- 

0.89 

Cover emmitance 

- 

0.85 

Flat plate absorber emmitance 

- 

0.31 

Cover reflectance 

- 

0.05 

Reflector reflactance 

- 

0.86 

Flate plate absorber reflactance 

- 

0.15 
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